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Abstract 

In this paper, we show theoretically that there exists quasi-liquid on the edges and vertices of 
snow crystals between — 4°C and —22° C, while the faces (0001) and (1010) have no quasi-liquid 
layers. Investigating the macroscopic theory of quasi-liquid and applying to the edges and vertices 
of the crystal, we find that the quasi-liquid becomes unstable above the critical supersaturation 
point, which is above the water saturation point. The thickness of this unstable quasi-liquid 
layer continues growing indefinitely. We interpret this behavior as corresponding to continuous 
production and overflow onto neighboring faces in a real system. We hypothesize that the unstable 
growth of snow crystals originates from the edges and vertices, and it is due to the overflow of 
quasi-liquid from the edges and vertices onto the neighboring faces, which are rough and lack 
quasi-liquid. Our hypothesis accounts for the qualitative behavior of the relations between the 
morphological instability and the water saturation in the snow phase diagram. 
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I. INTRODUCTION 



The snow crystal is one of the most beautiful things in nature, and it has attracted human 
interest since long ago. Its basic form is a hexagonal prism bounded by two basal (0001) 
and six prism faces (1010). Nakaya |I[ was the first to successfully produce snow crystals 
in the laboratory, and he investigated the relations between growth forms and experimental 
conditions (i.e., temperature and supersaturation relative to ice). Since his work, many 
experimental studies have been carried out 0. The snow phase diagrams obtained in these 
studies exhibit two important features (See Fig.|IJ). One is that three transitions occur in 
the basic crystal form: at temperatures near — 4°C (plates to columns), — 10°C (columns 
to plates), and — 22°C (plates to columns). These transitions correspond to changes in the 
growth rates of the faces (0001) and the faces (1010) with temperature |§. The second 
important point is that a morphological instability arises when the supersaturation relative 
to ice becomes sufficiently large. Especially above the water saturation point (saturation 
relative to supercooled water), the crystal growth originates from the edges and the vertices, 
and as a result, needles and sheaths are produced between — 4°C and — 10°C, and dendrites 
and sectors are produced between — 10°C and — 22°C 

The basic growth form of snow crystals is determined by the most slowly growing faces. 
The basal (0001) and prism (1010) faces grow slowly and become the bounding faces of the 
crystal. Assuming that the faces grow through the spreading motions of two-dimensional 
nuclei 0], and taking into account the vapor diffusion field surrounding the crystal, Frank 
proposed a possible explanation of the dendritic growth of snow crystals |J. Based on the 
same idea, Yokoyama and Kuroda || simulated the growth of an infinitely long hexagonal 
column (i.e. a two-dimensional snow crystal), and they succeeded in showing that six pri- 
mary branches are produced from the vertices of the column. The instability that results in 
the production of these branches is due to the nonuniformity of the vapor density surround- 
ing the polyhedron crystal. Their simulation, however, also showed that when secondary 
branches are produced, the size of the primary branches becomes comparable with that of 
the original hexagonal column. This is not a realistic modeling of the dendrites in actual 
snow crystals. Furthermore, this theory cannot account for the snow phase diagram dis- 
cussed above. It is thus apparent that we must take into account other, yet unknown effects 
to explain the dendritical instability. 
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FIG. 1: Variation of ice crystal form with temperature and supersaturation. Here Ap v is the 
vapor density excess. (Based on laboratory observations 0.) 

The three experimentally observed transitions in the basic crystal form at low supersat- 
uration were accounted by the theory of Kuroda and Lacmann |TT|. This theory is based 
on the well-known fact that the surface of ice existing just below 0°C is coated with a thin 
liquid like layer (quasi-liquid layer). This quasi-liquid layer has been observed using a va- 



riety of experimental techniques |L3|, [14], [15], [16], [L7|, |T§, [19| . Theoretically, Weyl || gave 
qualitative arguments for its existence, and Fletcher |J developed Weyl's arguments into a 
quantitative form. The starting point of Kuroda and Lacmann is a phenomenological model 
of the quasi-liquid layer proposed by Lacmann and Stranski [110]. On the basis of this model, 
Kuroda and Lacmann predicted that there exists quasi-liquid on the face (1010) at lower 
temperatures than on the face (0001). Recently, it has been observed in both experiments 



jig , |18|1 and molecular dynamics simulations |£0| that there is a difference in behaviour of 
the quasi-liquid layer on basal and prism faces. Kuroda and Lacmann hypothesized that the 
surface structure of ice changes with decreasing temperature from a surface covered with 
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a quasi-liquid layer to a rough surface without a quasi-liquid layer and finally to a smooth 
surface. According to their description, the changes in crystal form are due to the anisotropy 
in these surface structural transitions between basal and prism faces of ice. 

In this paper, we investigate the interaction between the quasi-liquid and the edges and 
vertices of snow crystals, whose growth originates from the edges and the vertices from — 4°C 
to — 22°C and above the water saturation point. In Section 2, we present a macroscopic 
model of a quasi-liquid layer and the explanation of the crystal form transitions proposed by 
Kuroda and Lacmann. In Section 3, we give an argument asserting that quasi-liquids remain 
on the edges and vertices in the temperature region between —4° C and — 22°C, where the 
faces (0001) and (1010) have no quasi-liquid layer. In Section 4, to examine curvature effects 
on the quasi-liquid, we investigate a model of an ice particle covered with a quasi-liquid layer 
in a vapor environment. In Section 5, we construct a macroscopic model describing a snow 
crystal whose six prism faces are without quasi-liquid but whose edges are covered with 
quasi-liquid. Using this model, we show that a novel instability arises on the edges and the 
vertices above the water saturation point due to the presence of the quasi-liquid. In Section 
6, we present a new description of the relation between morphological instability and water 
saturation in the snow phase diagram. In Section 7, we give a summary. 



II. QUASI-LIQUID LAYER ON A FLAT SURFACE 

A quasi-liquid layer can exist in a stable state because its presence reduces the surface 
free energy of the system. The energetic advantage of surface melting is represented by Act 
TOfl , defined as 



Aa = ct vs - a s£ - a iv , (1) 

where ct vs is the surface energy per unit area of a vapor-solid interface, a s e that of a solid- 
liquid one, and ai v that of a liquid-vapor one. This parameter is positive for a system of 
water and ice, and thus in this case the existence of the quasi-liquid lowers the surface energy. 



The free energy per unit area of a quasi-liquid layer of thickness 5 is given by M, 1 1 



AGpw^) = a £v + a s£ + AaW(S) + —(fi e - Ms), (2) 

where V q is the molecular volume in the quasi-liquid layer and fie and /i s are the chemical 
potentials per molecule for the bulk liquid and the bulk solid, respectively. Here the function 



W(8) must satisfy the conditions W(0) = 1 and W(oo) = 0, which result from the conditions 
AG p i ane (0) = a vs and AGpi ano (oo) = a ev + a s i + 5(p e - p s )/V q \s=oo- To this time, two types of 
W(S) have been used (see, for example, Ref.[pl|1). One is short range, Ws(<5) = exp(— S/A), 
and the other is long range, Wj^(5) = (1 + 5/A)~ n |TTJ, where n is a positive integer (n = 2 



for the Van der Waals forces [^T|]) and A is a parameter corresponding to the characteristic 
interaction length of the molecule in the quasi-liquid. Applying the minimization condition 
<9AGpiane/cM = 0, we find that the equilibrium thickness of the quasi-liquid layer is 

AQ m (T m — T) 



< = -An 



for W{8) = W s (8), and 



e q = -A + 



AaV q T m 



T 



(3) 



nA n AaV q ^ m ■ 

m \-L m ) j 



(4) 



for W{8) = Wl(5), where Q m is the energy of melting per molecule, T is the absolute 
temperature, and T m is the melting temperature. Here, we have used the equation fig — /i s = 
Qm(T m —T)/T m , which is derived from the Gibbs-Duhem relation assuming constant pressure 
and small T m — T. As seen in Eqs.@ and the thickness of the quasi-liquid layer decreases 
monotonically with falling temperature. The above equations should be considered valid only 
until the thickness becomes on the order of a few monolayers, but beyond this point, we 
can consider the layer to have vanished. At this temperature, the surface strongly adsorbs 
H 2 molecules, and it thus become rough at the molecular level. The number of adsorbed 
molecules decreases with falling temperature, until eventually it vanishes, and the surface is 
thus smooth. 

The parameter Act can be approximated in terms of the number density of the broken 
bonds per unit area p as 

where Q s is the energy of sublimation per molecule and Ab is the number of the bonds 
per molecule in the crystal p2| . For snow crystals, we know that pnoio) > P(oooi), from 
consideration of the surface molecular structures of ice, and Q s > Q m , from experimental 
tables [p^l - Thus A<7( 10 i ) > Act( 00 oi). According to Eqs.(|3]) and (f|), this implies that the 



quasi-liquid layer on a face (1010) remains at lower temperature than that on a face (0001). 



Thus snow crystals have the following four different surface structures [11|: (I) both the 



faces (1010) and (0001) are covered with a quasi-liquid layer; (II) the face (1010) is covered 
with a quasi-liquid layer, while the face (0001) has no quasi-liquid layer and is rough at the 
molecular level; (III) neither has a quasi-liquid layer, but (1010) is rough, while (0001) is 
smooth at the molecular level; (IV) neither has a quasi-liquid layer, and both are smooth at 
the molecular level. These structural differences result in differences in the growth speeds 
of these faces, since the mechanisms governing the growth depend on the surface structures. 
Assuming that the first surface structure transforms into the second at — 4°C, the second 
transforms into the third at — 10°C, and the third transforms into the fourth at — 22°C, 
Kuroda and Lacmann succeeded in describing the three crystal form transitions in the phase 



diagram [[LI] . 

The above consideration regards only the surface structures of faces. The situation be- 
comes more interesting when we also consider the surface structures of the edges and vertices. 
There exhibit transitions that differ from those of the basal and prism faces, as we shown 
in the following section. 



III. QUASI-LIQUID ON THE EDGES AND VERTICES OF SNOW CRYSTALS 

First, we approximate the energetic advantage represented by the surface melting pa- 
rameter for the edges and the vertices of snow crystals, A<T( e d g e & vertex) • In the process of 
snow crystal formation, the edges and vertices are formed as the intersections of metastable 
and unstable faces, such as (1120), (1011) and (1121). For this reason, we approximate 
Ac"(edge & vertex) by taking the average of Act for the intersecting faces. It is known from the 
molecular structure of ice that these metastable and unstable faces are rough and have large 
broken bond densities in comparison with the faces (0001) and (1010). Thus, from Eq.(|), 
we see that A<7( e d g e & vertex) is much larger than Acx(oooi) and A(j( W i y. 

Aa (oooi) < Aa ( 10 i ) < Act ( e d ge & vertex) • (6) 

This implies that for the edges and vertices, it is more energetically advantageous to have a 
surface transition from solid to liquid, and thus that they have a greater affinity for quasi- 
liquid. It is suggested by Eqs.(H) and flD that there exists quasi-liquid on the edges and 
vertices at lower temperatures than on the faces (0001) and (1010). 

In Fig|| we present surface structures of snow crystals in which the quasi-liquid of the 
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FIG. 2: The equilibrium thickness of a quasi-liquid layer for l^ s (<5) and W L (5) as functions of tem- 

o 

perature for various surface orientations. We chose A = 10 A, corresponding to the characteristic 



coherence length of a water molecule in the quasi- liquid |9|. This is also the cluster size in water [25]. 

o 

The quantity Act is computed using the monolayer thickness (3 A) at — 4°C , — 10°C and — 22°C. 
( From Eq.(|), we obtain Acr(OOOl) = 31 erg/cm 2 , Acr(1010) = 37 erg/cm 2 and Act (1120) = 54 
erg/cm 2 , where we use experimental values ae v = 76erg/cm 2 at 0°C, Q s = 8.48x 10~ 13 erg/molecule, 
and Q m = 1.0 x 10~ 13 erg/molecule J] lj . These values of Act are not accurate, because a vs for an 
actual relaxed vapor-solid interface is less than pQ s / '(2N&) for a freshly cut surface p2j| . Also, for 
actual values of a s £, we must take into account the fact that the water-ice interface has a diffuse 
structure throughout the thickness of several molecular layers [ p0[| . However, it is certain that 

A<7( 001) < A<7( 10 io) < A<7( e dge & vertex)-) 

edges and vertices are taken into account. We find the following: (I) from 0°C to — 4°C, all 
the faces, edges and vertices are covered with quasi-liquid; (II) from — 4°C to — 10°C, the 
faces (1010), edges and vertices are covered with quasi-liquid, while the faces (0001) have no 
quasi-liquid and are rough; (III) from — 10°C to — 22°C, the edges and vertices are covered 
with quasi-liquid, while the faces (0001) are smooth and the faces (1010) are rough, both 
without quasi-liquid; (IV) below — 22°C, no quasi-liquid exists. In the following sections, we 
examine the stability of the quasi-liquids on the edges and vertices in a vapor environment. 
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IV. QUASI-LIQUID LAYER ON A CURVED SOLID SURFACE 



To investigate the curvature effect on quasi-liquid at the edges and vertices, we consider 



here three- and two-dimensional ice particles covered with quasi-liquid [26]. We assume that 
the particles are in a vapor environment and have radius h + S, where h is the radius of 
the solid core and S is the thickness of the quasi-liquid layer. Although the main reason we 
consider this system is to examine curvature effects on a quasi-liquid, the system itself has 
meaning as a model of ice nucleation in a vapor environment. 
We assume a free energy of the above described system as 



AG particlc (M) =2 D -\ 



h D -\a sE + AaW{8)) + {h + 5) D - l a N 



~{(h + 5) D - h D } ^— ^ - }-h D ^-^ 

L) Vq L) V s 



(7) 



where D = 3 for three-dimensional particles and D = 2 for two-dimensional particles, fi v , 
fie and fi s are the chemical potentials per molecule for the bulk vapor, bulk liquid and bulk 
solid, and V q and are the volumes per molecule for the quasi-liquid and the solid. The 
first two terms on the right-hand side of Eq.(^|) are the surface energies of the system, and 
the last two terms are the bulk contributions. The surface energy contribution must coincide 
with 2 D - 1 nh D - 1 a vs for 5 = and 2 D - 1 nh D - 1 a si + 2 D - 1 n(h + 5) D - l o iv for 5 -> oo. Thus 
W(0) = 1 and W(oo) = 0. As in Eq.@, we use two types of W{5): W s {5) = exp(-5/A) and 
W L (5) = (1 + 5 1 A)- 2 . The quantity AG partidc (/i, 5)/(2 D - 1 7rh D - 1 ) reduces to Eq.@ when 
fi v = fi s and h — > oo. Since we are interested in the growth of ice in a vapor environment, 
we assume /i v > /i s vand fi£ > /i s in the following. 
The equilibrium condition is given by the equations 



<9AG par ticic 



0, (8) 



86 

0. (9) 



8 AG particle 



8h 

In equilibrium, we find a useful relation between 5 and h. From Eq.(|8]), we have 

(D -l)(h + 5) D - 2 a £v -(h + 6 )°-±Vlp^ = - h »-i Aa ^>l. (io) 

v q ad 

Thus, Eq.(||) becomes 

(D l)W(8) h dW ^ - _ ( t^-t^ V±zM P- 1 )^^ ril s 
1 ' 1 ' d5 ~ Aa V V q h )' 1 ' 
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This equation determines S as a function of h. If we use W(5) = Ws(5), we obtain 



S s (h) = -Aln 



h 



-Aln 



AX(h) 



(D-l)A + h_ 
and if we use W(5) = Wl(8), we obtain 

(D-l)Aa 1 



(12) 



6 L (h) = -A + A[ X (h,\(h)) + 



where X(h) and x(h, A) are defined by 



3hX(h) X (h,X(h)) 



(13) 



(14) 



and 



x(M) 



Act / / (£> - l) 3 A 2 Aa 
AX + 2WX 



(15) 



The above functions 5 s (h) and <5 L (/i) diverge at h = h c = (D — 
l)oW {(a* v — Hs)/V s — (a*v — and when > /i c , they monotonically decrease 

and approach the following 5 s (oo) and <5 L (oo) asymptotically: 

'AA(oo) 



<5 b (oo) 



5 L (oo) 



-Aln 



-A + 



Act 

2A 2 Acx 
A(oo) 



(16) 
(17) 



No real solution of 5 s (h) and £ L (/i) exists for < h < h c . 

When V s > V q , which holds for ice and water (and quasi-liquid), the following relations 
hold: 



h c >(D-l)V a <T gt /(n i -fi B ), 5 s (oo)>5 c s q , 5 L (oo)>5, 



'eq- 



(18) 



The quantities <5f q and 5^ q are given in Eqs.(|3]) and (f|). Thus we find the following two 
properties of h and S. One is that h > (D — l)K;Cr s ^/ (fie — A* s ) i n equilibrium. For D = 3, this 
means that the equilibrium radius of the solid part of the particle ,h, is always larger than 
the critical radius of ice nucleation in supercooled water. The other is that 5 s (h) in Eq.(|12"D 
and 5 L (h) in Eq.flT3|) are always greater than <5f q in Eq.(^|) and S^ q in Eq.(||), respectively. 
This implies that the equilibrium thickness 5 on a curved surface is greater than that on a 
plane surface, 6 eq . 



9 




supersaturation relative to ice ^(Pv)=(Pv-Ps)/Ps 

FIG. 3: The equilibrium states of a quasi-liquid on a curved surface as a function of the super- 
saturation relative to ice. For this figure, we set T = — 15°C, Act = 36.2 erg/cm 2 , A = 10 A and 

°3 

Vq = V s = 30A • From experimental values for water and ice, we have a s £ = 24.25 erg/cm and 



a iv = 78.25 erg/cm 2 at -15°C §| 



We solved Eqs.(|8]) and @ numerically. In Fig.|3|, we plot the ratio 8/(h + 5) in equilibrium 
as a function (P v — P s )/P s , where P v is the actual vapor pressure and Pi and P s are the 
equilibrium vapor pressures for the bulk liquid and the bulk solid, respectively. To obtain 
Fig|| we have used the relations 

/i v (T,P v ) = fx v {T) + kTlnP v , (19) 
IH(T, Pi) = fiy{T, Pi) = /i v (T) + kTlnPi, (20) 
/i s (T, P s ) = /i v (T, P s ) = /i v (T) + A;TlnP s , (21) 

where k is the Boltzmann constant and fJ> v (T) is the chemical potential of vapor at P v = 1. 
When Pi > P v > P s , the equilibrium state is uniquely determined. When P v > Pi, in 
addition to this state, another equilibrium state that has a thicker quasi-liquid layer appears. 
This state appears because above Pe it is possible not only for the solid but also for the liquid 
to grow in the vapor environment. As the supersaturation P v increases, these two equilibrium 
states approach each other, and at P v = P c they coincide. Above P c no equilibrium state 
exists. All the equilibrium states obtained here are labile equilibrium states. From the 
viewpoint of ice nucleation in a vapor environment, they all correspond to a critical radius 
of nucleation. 
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6 



FIG. 4: The quasi-liquid remaining on an edge. 

V. INSTABILITY OF THE QUASI-LIQUID REMAINING ON THE EDGE OF 
THE CRYSTAL 



Let us now consider a snow crystal whose edges are covered with quasi-liquid and whose 
faces have no quasi-liquid layer (see Fig.0). For the sake of simplicity, we consider only the 
six edges at which two different prism faces meet (see fig.f|). For other edges and vertices, 
similar analyses can be carried out and yield qualitatively similar results. 

The six prism faces (1010) grow slowly in comparison with the edges, and therefore we 
use the approximation that their speed of growth is zero. As shown in Fig.|4], we also simplify 
the shape of the quasi-liquid layer on the edges. This simplification is justified when 5 < h. 
For this system, the free energy of the snow crystal is given as 



AG prism (M)/12 



R 



°"(1010) 



R 2 ^ 



2^3 V B 



h 
IV3 



°"(1010) 



h 2 fi y 



2v/3 V s 



+ 



7T 

L6 



M°s<?(cdgc) + A<J( ed ge)W(5)) 



+ l~{h + S) + ^h 2 + (h + S) 2 - 2h(h + 8) cos 



~{ T ^h + 5) 2 -h 2 ) 



7T 



-—k 
12 



,(22) 



where cx(ioio) is the surface energy for a prism face that has no quasi-liquid layer and adsorbs 
H 2 molecules, and cr S £( e d g c) and A<j( e d ge ) are the surface energy of the solid-liquid interface 
and the energetic advantage of surface melting (given by Eq.(||)) for the surface of the edge, 
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respectively. The curvature effects considered in the previous section are taken into account 
by the terms inside the last set of square brackets on the right-hand side of Eq.fl22|). When 
h > 5 and (£) 3 ~ 0, Eq.(^) becomes 



AG pri sm(M)/12 



where 



+ 



R R 2 /x v — jj, s 

7T 

g^(0"s«(edge) + A(T( ed g e )W / ((5)) 



h h 2 fi v — fM s 



9 V 2| 2h 4 1 V36 



h 



h } a lv 



36 V V n 12 K 



(23) 



h 2 + (h + 5) 2 - 2h(h + 5) cos(^) = J5 2 + 4h(h + 5) sin 2 (^) 



2sm(^{l + i- + l( sin 



-2, 



, 7T , 

36' 



and we have used the approximations sin (y|) ~ ^ and sin (JU ~ ^. The equilibrium 
condition is given by 

dAG r 



' prism 



88 



0. 



8AGprism p. 



(24) 
(25) 



From Eq. (fZ5j) , we find 



i 



73° (1010) - 6 U vs(edge) 

n — n — — 7 r 

V3 el y s 



(26) 



where o" vs ( e d g e) is the surface energy of the vapor-solid interface for the edge. In order 
for foo to be positive (which means that the curved edge is energetically stable) and for 



( dAG v 



h=ho 



< to hold (which means that a curved edge with a quasi-liquid layer is 



V 96 J 5=0 
energetically stable), cr(ioio) must satisfy 



a 



(1010) 



> 



v/3 



7T 



~C"vs(edge) 4' 



5 / _1 7T \ /Jy -^s 

V3 6 y °^ y s 
An-, , . 9W(0) , 5 Mv~M£ 

-^o( edge ) 9J -r g y q 



(27) 



This equation gives a lower bound on a 



(1010)- 
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FIG. 5: The stability of the quasi-liquid remaining on an edge as a function of the supersaturation 
relative to ice. We plot the case of W s (5). For W^(5) , qualitatively similar results are obtained. 
The vertical solid line at supersaturation 0.167 indicates the water saturation point. The horizontal 
line at Sf q indicates the equilibrium thickness of a quasi-liquid layer on a planar surface, given by 
Eq.(@). We used T= -15°C (258.15K), cr (10i0) = 129.82erg/cm 2 and Acj (edge) = 36.2erg/cm 2 . We 
also used the experimental values at — 15°C, <7 s f( e dge) = 24.25erg/cm 2 and C£ v = 78.25erg/cm 2 
p2fl (and therefore 0" vs ( e d ge ) = 138.7erg/cm 2 ). This model is justified when S is larger than the 
mono-layer thickness and h is sufficiently large. Thus we ignore the unstable solutions existing for 
5 < 3 A and h < 30 A. This type of figure can be obtained when (T(ioio) > 128.9erg/cm . When 
/i v — /jl s — ► oo, the right-hand side of Eq,(^7|) converges to 133.07 erg/cm 3 . This value is reasonable 
for a rough face (1010) having no quasi-liquid layer, (c.f. c vs (ioio) = pQs/(2A r B) = 128erg/cm 2 , 
where Q s = 8.5 x 10 _13 erg/molecule ). The critical supersaturation C(-fc) is 0.265. 
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Numerical solutions of Eqs.fl24|) and fl25|) are plotted in Fig.|5]. (Here we have used 
Eqs.([T9|)-(|21"|) again.) Remarkably, we find a stable equilibrium state when P v < P c . In 
this case, a snow crystal with small 5 and h evolves toward a stable state, at which point it 
ceases growing. 

When P v > P c (> Pi), no stable solution satisfying Eqs.fl2"4]) and (^Sj) is found. It is 
interesting that the thickness 5 continues to grow in this case. This suggests that a novel 
instability arises on the edges of snow crystals when P v > P c . Of course, in a real system, 5 
cannot continue to grow indefinitely, and thus this result reveals a limitation of our model. 
Our model breaks down when S becomes comparable with h, and our simplification regarding 
the shapes of the quasi-liquids on the edges is no longer justified. In a real system, as the 
thickness of the quasi-liquid grows, eventually it begins overflowing onto the neighboring 
faces (1010) or (0001) to reduce the surface energy. We interpret the continuous growth of 6 
in our model as corresponding to the continuous overflowing of quasi-liquid in a real system. 
If the neighboring face is rough at the molecular level, the overflowing quasi-liquid sticks 
to the neighboring face, and immediately turns to be solid near the edge. As we show in 
the next section, this interpretation allows our model to explain the phase diagram of snow 
crystals. 

VI. PHASE DIAGRAM 

In this section we give a new description of the snow phase diagram. Let us reconsider 
the discussion of Sectionp|. In the temperature regions (I), 0°C to — 4°C, and (IV), below 
— 22° C, the quasi-liquid on the edges and vertices does not play a special role. Thus we 
focus on the temperature regions (II) and (III), — 4°C to —22°. 

In the temperature region (II), — 4°C to — 10°C, the six prism faces (lOlO), all the edges 
and all the vertices are covered with quasi-liquid, while the two basal faces (0001) are rough 
at the molecular level and have no quasi-liquid layers. As seen in Fig. |5], when the vapor 
pressure P v is low, the quasi-liquids on the edges and vertices are stable. In this case, the 
edges and the vertices do not play a special role in the snow crystal formation. However, 
when P v is higher than a certain value P c (> Pi), the quasi-liquid on the edges and vertices 
overflows onto the neighboring faces. The quasi-liquid overflowing onto (lOlO) spreads over 
the quasi-liquid layer, while the quasi-liquid overflowing onto rough (0001) sticks near the 



14 



edges and vertices and turns to solid. In this case, two types of overflow are possible. One 
is overflow from the vertices, and the other is overflow from the edges where the basal and 
prism faces meet. If the former becomes the main type of overflow, the resultant snow 
crystal should grow like a needle. If the latter becomes the main overflow according to our 
model, the resultant snow crystal should grow like a sheath. 

In the temperature region (III), — 10°C to — 22°C, the edges and the vertices are covered 
with quasi-liquid, while the basal (0001) and prism (1010) faces possess no quasi-liquid layer. 
The basal faces are smooth and the prism faces are rough at the molecular level. Here again, 
as seen in Fig. |, when P v is low, the quasi-liquid is stable, and thus the snow crystal grows 
according to the theory of Kuroda and Lacmann. When P v is higher than P c (> Pi), the 
quasi-liquid overflows onto the neighboring faces. The quasi-liquid overflowing on to smooth 
(0001) spreads over (0001) and turns to solid, while the quasi-liquid overflowing onto rough 
(1010) sticks near the edges and vertices and turns to solid. In this case, three types of 
overflow are possible. The first is overflow from the vertices, the second is overflow from the 
edges where the basal and the prism faces meet, and the third is overflow from the edges 
where two prism faces meet. If the amounts of quasi-liquid involved in all overflows are of 
the same order, the resultant snow crystal should grow like a sector at about — 10°C. If the 
first and (or) the third overflow are the main types of overflow, the resultant snow crystal 
should grow like a dendrite at about — 15°C. Contrastingly if the first and (or) the third 
overflow are the main overflows, but the surface structure of prism face changes from a rough 
surface to a smooth surface with temperature falling, the resultant snow crystal should grow 
like a sector again. This is because the continuous overflow spreads over smooth prism faces 
and turns to solid successively. In this case, the overflow can act as a step source for prism 
faces and thus the snow crystal grows in the manner hypothesized by Frank || and found 
in the simulation of Yokoyama and Kuroda || . We believe that this type of sector can be 
observed at about — 22° C. 

We note that the critical supersaturation shown in Fig.[5]is 0.265 at — 15°C, which agrees 
with the experimental supersaturation, for which dendrites are observed above 0.20 (see 
Fig.©. 
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VII. SUMMARY 



In this paper, we have shown theoretically that quasi-liquid layers exist on the edges 
and vertices of snow crystals between — 4°C and — 22°C, for which temperatures there is no 
quasi-liquid on the faces (0001) and (lOlO). 

To investigate the curvature effect on quasi-liquid layers at the edges and vertices, we 
considered an ice particle coated with a quasi-liquid layer, and derived the equilibrium over 
a range of supersaturation values. We found that below the water saturation point, the 
equilibrium state is uniquely determined. However, above this point, another equilibrium 
state exists because it is possible for bulk supercooled water to exist. Above the critical 
supersaturation point, no equilibrium state exists. We showed that the quasi-liquid layers 
for equilibrium states of a particle are thicker than the quasi-liquid layer for the equilibrium 
state of a planar surface. 

Next, we considered a two-dimensional crystal. We found that there are temperatures for 
which the six prism faces possess no quasi-liquid, while the edges and vertices are covered 
with quasi-liquid. We have examined the stability of this system. We found that below the 
critical saturation, which is higher than the water saturation, the quasi-liquid layer exists 
in a stable state on the edges. However, above the critical supersaturation, the quasi-liquid 
becomes unstable, and continues to grow indefinitely. We interpret this indefinite growth as 
implying overflow onto the neighboring faces in a real system. 

On the basis of our results, we have proposed a new description of snow crystal growth 
according to which the unstable growth of snow crystals, which prefer the edges and the 
vertices, is due to the overflow of the quasi-liquid from the edges and the vertices onto 
neighboring faces that are rough and have no quasi-liquid layers. We have shown that this 
overflowing occurs above the water saturation point and that these surface conditions of 
snow crystals are realized between — 4°C and — 22°C. This description naturally accounts 
the relation between the morphological instability and the water saturation in the snow 
phase diagram. 

While we have succeeded in explaining the snow phase diagram qualitatively, for a 
more complete understanding, the details of crystal growth should be clarified by three- 
dimensional simulation. We also would like to show that secondary branches of the den- 
drites are produced by the overflow of quasi-liquid from the edges and vertices of the primary 
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branches. These are future problems. 

In this paper, we have focused on the quasi-liquid of snow crystals. However surface 
melting is observed in many classes of solids, including metals, semiconductors, solid rare 
gases [BTJ. The model and method presented in this paper are also applicable to these 



materials, and should be useful for the investigation of general properties of surface melting. 
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